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ABSTRACT 

When  electron-electron  correlations  are  important,  it  is  often  necessary  to  use  “exact”  numerical 
methods,  such  as  Lanczos  diagonalization,  to  study  tl  o  full  many-body  Hamiltonian.  Unfortunately, 
such  exact  diagonalization  methods  are  restricted  to  small  system  sizes.  We  show  that  if  the  Hubbard 
U  term  is  replaced  by  a  “periodic  Hubbard”  term,  the  full  many  body  Hamiltonian  may  be  exactly 
solved,  even  for  very  large  systems,  though  for  low  fillings.  However,  for  half-filled  systems  and  large 
U  this  approach  is  not  only  no  longer  exact,  it  no  longer  improves  extrapolation  to  larger  systems. 
We  discuss  how  generalized  “randomized  variable  averaging”  (RVA)  or  “phase  randomization”  schemes 
can  be  reliably  employed  to  improve  extrapolation  to  large  system  sizes  in  this  regime.  This  general 
approach  can  be  combined  with  any  many-body  method  and  is  thus  of  broad  interest  and  applicability. 

INTRODUCTION 

Solutions  of  the  Peierls-Hubbard  Hamiltonian  (PHH)  [1]  for  small  clusters  are  strongly  dependent 
on  the  boundary  condition  (BC);  i.e.,  for  the  case  of  polyacetylene,  whether  the  carbon  atoms  are 
viewed  as  being  on  a  real  ring,  such  as  benzene,  a  real  chain,  such  as  (l,3,5)-hexatriene,  or  a  more 
exotic  geometry,  such  as  anti-periodic  BCs.  We  have  found  that  [2,3]  to  reduce  finite  size  corrections 
to  the  calculated  optical  gaps  and  spectra,  thus  improving  the  extrapolation  to  the  infinite,  case,  it  is 
effective  to  employ  novel  “averaging”  techniques  to  “randomize”  the  many-body  levels  which  influence 
those  correlation  functions  of  interest.  Briefly,  by  randomization,  we  mean  simply  that  the  total  energy 
of  the  system  is  to  be  viewed  as  a  weighted  average  of  the  energies  derived  separately  for  each  of  several 
different  “random”  values  (RV)  of  a  particular  parameter,  which  may  be  the  BC,  a  hopping  integral, 
an  on-site  energy,  or  the  Hubbard  V:  E  —  XRV E{RV],  where  is  a  normalized  weighting 

factor  with  Y2{rv}  xRV  =  1.  Similarly,  any  correlation  function  can  be  viewed  as  the  average  of  of 
its  values  from  the  separate  RVs.  The  number  of  RV,  Nr y,  may  be  varied,  and  the  choice  of  the  RV 
and  xrv  may  be  truly  random,  or  may  be  a  pre-specified  set  of  values,  which  may  in  turn  depend  on 
the  quantity  being  studied.  The  case  of  phase  randomization  and  its  relation  to  Bloch’s  theorem  are 
described  below.  More  generally,  this  randomized-variable-averaging  (RVA)  technique  may  be  viewed 
as  a  procedure  to  mimic  disorder. 


We  focus  iu  this  manuscript  on  the  calculation  via  exa<  t  l.anczos  diagonajizat  ion  of  the-  : « -i I 
consistent  uniform  dimerization  and  one  photon  optical  absorption  of  the  ground  state  within  the  1  j). 
l/2-filled,  one-band  PHIL  However,  RVA  is  equally  applicable  to  the  calculation  of,  <  .g.,  the  phonon 
inodes  of  the  ground  state,  as  well  as  the  geometry  and  self-consistent  absorptions  of  doped  or  hiphei 
lying  (e.g.,  triplet)  states  [3],  luminescence  spectra,  structure  factors,  etc.,  in  single-  or  multi  band 
models  in  1-,  2-,  and  3-D,  via  Lanczos  diagonalization,  as  well  as  other  numerical  procedures  on  small 
lattices,  such  as  Monte  Carlo. 


PHASE  BOUNDARY  CONDITION  AVERAGING  AND  BLOCH’S  THEOREM 

Bloch’s  theorem  tells  us  that  the  single  particle  wavefunctions  of  a  system  of  size  M-N  with 
periodicity  N  and  periodic  boundary  conditions  (PBCs)  are  exactly  the  single  particle  wavelunctions 
of  the  ensemble  of  systems  of  size  N  for  each  of  the  phase  BCs  tp(N  +  1)  =  e‘^>(l),  <p  -  2nC/M . 
£  —  I,...,  M.  Restated,  the  (exact)  properties  of  a  larger  (single  parti'0'')  system  may  be  found  in¬ 
forming  a  (particular)  average  over  smaller  systems  with  different  BCs.  This  “boundary  condition 
averaging”  (BCA)  or  “phase  randomization”  is  a  special  case  of  our  general  RVA  procedure.  Blochs 
theorem  may  be  generalized  to  many  particle  wavefunctions  and  is  applicable  here  if  the  form  of  the 
Hamiltonian  is  modified  slightly.  The  one-dimensional  PHH  is 

H  —  ^(-*0  +  «^r)(c]  aCl+]  a  +  c\+l  a)  +  U^Tneintl  -f  -  A'  Sj  (1) 
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then,  if  the  original  problem  had  periodicity  N ,  («*,<,}  =  (ne+l, -N,a),  we  have,  in  a  mean-field  sense, 
the  “same”  Hamiltonian.  With  this  new  “periodic-Hubbard”  Hamiltonian  (periodic  PHH)  the  Bloch 
analysis  of  a  given  large  (periodic)  system  (N  ■  M  sites,  Ne  ■  M  electrons)  is  accomplished  by  considering 
a  small  system  ( N  sites,  Ne  electrons)  with  several  different  BCs.  The  method  scales  linearly  with  A/, 
allowing  one  to  handle  reasonably  large  systems,  and  sets  a  sound  theoretical  basis  for  the  empirical 
observation  (2]  that  RVA  can  be  used  to  smooth  optical  absorption  spectra  obtained  via  exact  finite-size 
diagonalization. 

Using  the  periodic-Hubbard  term,  Eq.  (2),  assuming  the  lattice  distortion  has  the  same  periodicity. 
£>{£  +  N)  =  S(i),  and  using  £^1,  =  Nc-  M,  one  can  show: 
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Note  that  H(iyj,k,l;$i,Sj)  is  independent  of  /<,  r\  and  M .  The  Bloch  analysis  on  the  many  body 
eigenfunctions  using  the  symmetries  of  H (i,j,  k,l;  * i,5j)  leds  to  a  4’  of  the  form: 


N 

=  E  E 

cr,  t,i  ni . nF  =  l 

M  M 

X  C"i+ni  N,<r,e  X  X  E  C"p  +  ep  N,ffpe 

<<i  =  l 


(1) 


P  ^P 


where  4>kiy...,kt . ,  «P)  =  . *,("1 ,  -,«<+/**•  np) ,  fc,  =  (2 vj/M),  je{0, ...,  A/ 

1},  and  p  =  Ne  ■  M.  We  find  for 
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(vq  comes  into  play  for  operators  like  c^N  ac^+\<3).  From  Eq.  (5)  we  see  Hper  is  diagonal  in  the  fc’s, 

Hver  =  Ep  ® kj,...,kp  >  (C) 


as  are  p  =  £/<r  »/,a  and  J  =  t  £*><r(-t0  +  a<5*)(c] i(7Q+t,<r  -  cJ+1(Tc<r,a),  and  the  eigenfunctions  of  7/pe.r 
can  be  written  in  the  form  Eq.  (4). 

The  symmetry  of  ^  implies  that  if  atl  -  at7  and  -  kt7  then  n(l  /  n(7 .  Thus  there  are  at  most 
2 N  of  the  k(  the  same.  We  postulate  that  the  ground  state  lies  in  the  manifold  with  each  of  the  M 
distinct  values  for  kt  occurring  Ne  times  and  (for  the  half-filled  band)  with  equal  numbers  of  up  and 
down  spins  (in  general,  the  largest  manifold).  This  can  be  checked,  and  we  stress  that  the  analysis  up 
to  this  point  is  exact  for  the  periodic  PHH.  In  this  manifold  we  can  write  as: 
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where  qt  =  (2 x£/M)  and 
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where  TVj  =  JVj  =  [[7£l]]  ([[  ]]  denotes  integer  value  of).  Note 
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and  so  solving  the  /7e  electron  problem  on  the  full  N  -M  sites  with  PBCs  for  fixed  q  is  exactly  equivalent 
to  solving 


per  'r  q  ,cr 


Hve  r  9 


PN*  \hN* 
q,Ce  q,a 


(9) 


on  N  sites  with  the  ^-dependent  BC  defined  by  Eq.  (8).  Thus,  to  find  the  exact  Ne  electron  eigenfunc¬ 
tions  on  N  ■  M  sites,  we  only  need  to  solve  Eq.  (9)  with  N  sites  and  NF  electrons  for  each  of  the 

KA  R  c  rtif  Kor  tKm  f  ho  full  Af  .  A/f  cit/a  rv  r/^Klorvi 


To  obtain  the  exact  eigenfunctions  of  the  full  Ne-M  electron  problem,  w«-  need  to  solve  for 
4>^(ai , a/if ).  This  is  where  an  approximation  must  be  made  to  be  able  to  solve  the  problem  numer¬ 
ically,  as  we  assume  N  is  already  as  large  as  computationally  feasible,  i  he  full  eigenvalue  problem 


yields 


To  solve  Eq.  (6)  we  must  use  some  approximation  for  Eq.  (10),  such  as  perturbation  theory  or  mean- 
field.  Thus  BCA  means  we  treat  electron-electron  correlations  within  the  A',  electron  manifold  exactly, 
and  electron-electron  correlations  between  Ne  electron  manifolds  approximately.  The  BCA  results 
reported  here  are  zero-order  perturbation  theory:  we  have  assumed  4>p(ai,  to  be  a  product  of 

^-functions;  i.e.,  for  the  ground  state  we  use: 
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Note  that  if  we  had  been  interested  in  systems  at  low  filling  ( Ne  electrons),  rather  than  near  half  filling 
(Ne  ■  M  electrons),  there  would  be  no  need  to  make  any  approximations. 


Comparison  of  the  periodic  and  standard  Peierls-Hubbard  Models  at  half-filling 

At  U— 0  Bloch’s  theorem  for  the  many  particle  wavefunctions  as  formulated  above  is  exact  at  all 
fillings.  For  the  half-filled  PHH,  we  have  just  shown  that  phase  BCA  involves  two  approximations: 
first,  replacing  the  Hubbard  term  by  Eq.  (2),  and,  second,  ignoring  correlations  between  the  Ne  particle 
wavefunction  leading  to  Eq.  (11).  We  now  test  these  approximations.  In  Fig.  1  we  show  the  minimum 
energy  dimerization  of  an  8-site  system  as  a  function  of  U  for  (i)  phase  averaged  solution  and  (it) 
the  exact  solution  of  the  half- filled  band  with  the  periodic  Hubbard  term  and  periodicity  2,4,  and 
8  (periodicity  8  is  the  the  usual  Hubbard).  The  approximate  periodic  (BCA)  solution  lies  between 
the  exact  periodic  solution  and  the  exact  pure-Hubbard  solution.  Since  the  original  problem  was  for 
U  •  nqnq  on  N  ■  M  sites,  one  might  argue  that  the  approximate  phase-averaged  result  for 

U  2_,<=1  on  N  sites  may  be  closer  to  the  desired  answer  than  the  exact  periodic-Hubbard  N  -  M 

site  result. 

In  Fig.  2  we  show  the  dimerization  amplitude  obtained  by  phase  BCA.  We  see  that  the  extrapolated 
infinite  behavior  at  small  to  intermediate  U  is  well  approximated  after  phase  averaging  for  Ar  =  10  even 
without  extrapolation.  The  agreement  is  best  for  small  {/,  while  no  change  from  the  pure  periodic 
behavior  is  found  for  large  U .  That  this  particular  technique  is  expected  to  have  no  effect  at  large  U 
can  be  seen  by  examining  the  effective  spin-Peierls  Hamiltonian  [4],  which  is  asymptotically  independent 
of  the  phase  of  the  BC.  However,  the  phase- averaged  analysis  does  lead  to  the  interpretation  of  the 
finite-size  result  as  the  result  for  the  larger  system  with  an  effective  periodic-Hubbard  interaction. 
Thus,  if  U  is  scaled  by  the  system  size,  as  is  done  in  Fig.  2,  the  infinite  behavior  should  be  more  easily 
extrapolated.  It  appears  that  this  is  at  least  approximately  true. 

To  show  how  phase  BCA  affects  the  optical  spectrum,  in  Fig.  3a.  we  show  the  (Lorentzian  broad¬ 
ened)  spectrum  at  U jt^—QA  for  N - 8  and  5  phase  BCs  (enough  to  yield  a  smooth  spectrum  at  (7=0). 
It  is  clear  how  the  different  BCs  “fill  in”  the  spectrum.  Fig.  3b  shows  the  spectrum  at  f//t0=;1.6.  Here, 
despite  the  several  BCs  used,  the  spectrum  remains  sparse,  due  to  the  loss  of  effectiveness  of  the  phase 
BCA  scheme  as  discussed  above.  If  we  were  to  use  a  broad  enough  Lorentzian  to  smooth  this  spectrum 


8U/N 


Fig.  1  (left).  Comparison  of  the  approximate  and  exact  periodic-Hubbard  minimum  energy  dimerization 
as  a  function  of  U  for  <0=2.5  eV,  a=4.1  eV/A,  and  A=21  eV/A2. 

Fig.  2  (right).  Dimerization  vs.  U  for  various  N:  phase  BCA.  Parameters  as  Fig.  1. 


Fig.  3.  The  Lorentzian  broadened  optical  absorption  spectra  calculated  using  the  phase  BCA  method 
for  (a)  U=l  eV  and  (b)  U- 4  eV.  Note  the  “sparseness”  of  the  spectrum  in  the  latter  case  indicating 
the  failure  of  the  phase  BCA  approach.  For  comparison,  we  also  show  the  spectrum  at  V  -  10  eV  using 
(c)  amplitude  BCA  and  (d)  an  RVA  procedure  where  U  is  varied  on  one  site.  Note  these  RVA  spectra 
are  dense  and  show  the  same  features.  Parameters  as  Fig.  1  but  <5=0.14  A. 


AMPLITUDE  BOUNDARY  CONDITION  AVERAGING 

The  modification  of  the  phase  of  the  BC  discussed  above  can  also  be  viewed  as  passing  a  random 
magnetic  flux  through  an  ensemble  of  closed  ring  and  studying  average  properties.  In  this  sense  we  have 
“randomized”  the  locations  of  the  momentum  space  states.  We  have  already  indicated  that  RVA  can 
help  us  transcend  the  limitations  of  this  “bond  phase/magnetic  flux”  approach.  We  can  “randomize” 
electronic  properties  by,  e.g.,  changing  a  local  hopping  or  an  on-site  energy  or  a  Coulomb  repulsion 
somewhere  on  the  chain.  One  could  also  introduce  an  additional  field  and  vary  it  about  zero.  Although 


[5],  there  is  as  yet  no  provably  accurate  prescription  for  for  arbitrary  E  ami  1.  However.  <miatn 
intuitive  rules  must  guide  us.  First,  whatever  change  is  made  to  randomize  must,  of  course,  do  so 
effectively.  Second,  the  change  in  the  system  must  be  negligible  as  the  lattice  size  is  increased  to 
infinity;  for  example,  if  only  one  bond  or  site  is  varied  from  calculation  to  calculation,  then  the  effect 
of  such  a  change  is  immaterial  in  the  thermodynamic  limit.  Finally,  the  behavior  for  small  lattice 
sizes  must  be  illustrative  of  the  infinite-size  limit.  Put  another  way,  one  must  still  be  able  to  make  ;i 
reasonable  extrapolation  to  the  infinite  chain. 

We  have  found  [2]  an  “amplitude  BCA”  (“scaled-hopping”)  technique  to  be  effective  and  pro¬ 
duce  results  in  good  agreement  with  expectations  based  on  both  strong-  arid  weak- coupling  analytic 
arguments.  For  this  method,  we  randomize  by  varying  the  magnitude  rather  than  the  phase  of  the 
“boundary”  hopping  -  i.e.,  between  sites  1  and  N  -  typically  from  -t0  to  -Mo  in  ten  to  twenty  equal 
steps.  The  individual  spectra  are  then  added  together  with  weights  x,  chosen  to  minimize  the  total 
length  of  the  final  curve,  though  giving  each  spectra  equal  weight  yields  virtually  identical  results.  This 
clearly  incorporates  the  special  case  of  using  only  periodic  (ar=l)  or  antiperiodic  (x--l)  rings,  the 
JT/nJT  difference  being  important  for  weak  coupling.  It  also  incorporates  the  case  of  the  open  chain 
(z=0),  which  has  the  “best”  single  BC  size  dependence  (though  with  attendant  “end  effect”  problems). 
Finally,  it  works  in  the  strong- coupling  limit. 

In  Fig.  3c,  we  show  the  spectra  produced  by  the  amplitude  BCA  technique  for  for  a  larger  Li  than 
that  in  Fig.  3b  where  phase  BCA  failed.  In  Ref.  [2]  we  showed  that  not  only  are  the  gross  features 
of  the  spectra  obtained  by  amplitude  BCA  in  agreement  with  strong-coupling  calculations  [6],  but  in 
addition  they  show  substantially  more  interesting  detail,  such  as  the  “decoupled- dimer”  peak  located 
at  ~  (7/2  +  y/(U/2f  +  (2t0)2  in  systems  with  strong  electron-phonon  and  electron-electron  couplings. 
To  emphasize  the  generality,  in  Fig.  3d  we  show  the  spectrum  obtained  by  varying  U  on  the  first  site 
between  0  and  20  eV.  This  large  U  variation  produces  states  in  the  optical  gap,  but  otherwise  the 
spectrum  is  unchanged  from  Fig.  3c.  We  feel  that  these  results,  and  those  in  Ref.s  [2,3],  demonstrate 
that  RVA  is  an  effective  means  for  reducing  (and  in  some  cases  practically  eliminating)  finite  size 
dependence,  yielding  results  that  can  be  confidently  extrapolated  to  the  infinite  size  limit,  and  which 
are  in  good  agreement  with  known  analytic  results. 

Acknowledgements.  We  wish  to  thank  J.  Bronzan,  J.  Gubernatis,  S.  Mazumdar,  and  D.  Scalapino, 
among  others,  for  many  useful  discussions.  JTG  was  supported  by  a  NRC-NRaD  Research  Associate- 
ship  through  a  grant  from  ONR.  Computational  support  was  provided  by  LANL  and  TMC. 

REFERENCES 

1  For  a  survey  of  the  PIIH  applied  to  conducting  polymers  see  D.  Baeriswyl,  D.K.  Campbell,  and 
S.  Mazumdar,  in  H.  Kiess  (ed.),  Conducting  Polymers ,  (Springer,  1991). 

2  E.Y.  Loh,  Jr.,  and  D.K.  Campbell,  Synth.  Metals,  27  (1988)  A499;  D.K.  Campbell.  J.T.  Gammol, 
and  E.Y.  Loh,  Jr.,  Int.  J.  Mod.  Phys.  D ,  3  (1989)  2131;  E.Y.  Loh,  Jr.,  D.K.  Campbell,  and 
J.T.  Gammel,  in  D.  Baeriswyl  and  D.K.  Campbell  (ed  ),  Proc.  of  the  NATO  ARW  on  Interacting 
Electrons  in  Reduced  Dimensions,  Torino,  Italy,  Oct.  3-7,  1988,  (Plenum,  1990). 

3  D.K.  Campbell,  J.T.  Gammel,  H.Q.  Lin,  and  E.Y.  Loh,  Jr,  Synth.  Metals,  49,  (1992),  in  press; 
D.K.  Campbell,  J.T.  Gammel,  and  Il.Q.  Lin,  this  conference. 

4  J.E.  Hirsch,  Phys.  Rev.  Lett.,  51  (1983)  296. 

5  R.  Julien  and  R.M.  Martin,  Phys.  Rev.  B,  26  (1982)  6173;  A.M.  Oles;  G.  Treglia.  D.  Spanjaard. 
and  R.  Julien,  Phys.  Rev.  B,  32  (1985)  2167. 

6  S.K.  Lyo  and  j.-P.  Gailinar,  J.  Phys.  C,  10  (1977)  1693;  S  K.  Lyo,  Phys.  Rev.  B.  18  (1978)  1851; 
J.-P.  Gailinar,  J.  Phys.  C,  12  (1979)  L335. 


